arXiv:hep-lat/0211020v2 13 Feb 2003 


CERN-TH /2002-304 
CPT-2002/P.4442 
DESY-02-178 
IFIC/02-52 
FTUV-02-1111 
hep-lat/0211020 


Finite-Size Scaling of Vector and Axial Current Correlators 

P.H. DamgaarcP, P. Hernandez 13 ’ 1 , K. Jansen c , M. Laine b , L. Lellouch d 


a Niels Bohr Institute, Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 
b Theory Division, CERN, CH-1211 Geneva 23, Switzerland 
c NIC/DESY Zeuthen, Platanenallee 6, D-15738 Zeuthen, Germany 
d Centre Physique Theorique, CNRS, Case 907, Luminy, F-13288 Marseille, France 


Abstract 

Using quenched chiral perturbation theory, we compute the long-distance behaviour of two- 
point functions of flavour non-singlet axial and vector currents in a finite volume, for small 
quark masses, and at a fixed gauge-field topology. We also present the corresponding pre¬ 
dictions for the unquenched theory at fixed topology. These results can in principle be used 
to measure the low-energy constants of the chiral Lagrangian, from lattice simulations in 
volumes much smaller than one pion Compton wavelength. We show that quenching has a 
dramatic effect on the vector correlator, which is argued to vanish to all orders, while the 
axial correlator appears to be a robust observable only moderately sensitive to quenching. 
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1. Introduction 


The low-energy dynamics of QCD is governed by a spontaneously broken chiral symmetry. 
Thus, close to the chiral limit, the correlation length associated with the Goldstone bosons 
is very large. In a finite volume, as required by lattice simulations, this correlation length 
can easily overtake the linear extent of the box. When this happens, the zero-momentum 
modes of the Goldstone bosons can no longer be treated perturbatively and finite-size effects 
become important. 

The general procedure for describing this situation in chiral perturbation theory is the 
e-expansion of Gasser and Leutwyler |T (see also [2j). It yields precise predictions for the 
volume and quark-mass dependences of long-distance observables in terms of a few infinite 
volume low-energy constants. The comparison of these predictions with the volume and mass 
dependences of the same observables computed in lattice QCD then permits the extraction 
of the constants. Such lattice studies have become possible recently thanks to the advances 
with Ginsparg-Wilson formulations [3j of lattice fermions which implement a continuum-like 
chiral symmetry at finite lattice spacing 

In full QCD, all the relevant two-point functions of currents and densities have already 
been calculated in the e-regime UU- However, due to the numerical cost of Ginsparg- 
Wilson fermions, lattice simulations in this regime are currently restricted to the quenched 
approximation. It is therefore important to evaluate the modifications brought about by this 
approximation. In a previous paper scalar and pseudoscalar correlators were calculated 
in the quenched approximation, in sectors of fixed topology. We extend that work here by 
presenting the calculation of flavour non-singlet axial and vector correlators. 

At leading order in the chiral expansion, the chiral Lagrangian is parametrized by two low- 
energy constants, the chiral condensate X and the pion decay constant F. The value of the 
first one can be accurately extracted from the spectral density of the low-lying eigenvalues 
of the Dirac operator (see, e.g., ref. (SD- It has also been obtained by measuring directly the 
quenched finite-volume quark condensate m-m- To determine F, on the other hand, we note 
that because it quantifies the strength of the coupling of Goldstone boson fluctuations to the 
vacuum, it is natural to consider correlation functions of these fluctuations. As shown in U El, 
the scalar and pseudo-scalar correlation functions are proportional to X 2 at leading order, 
with F appearing at the next-to-leading order. However, in the quenched approximation 
these correlators also depend on additional unphysical constants (m,Q and a), associated with 
the flavour singlet field [Zj, which clearly makes numerical determinations rather difficult, 
while maybe not impossible m- 

By contrast, the axial and vector current correlators are proportional to F 2 at leading order, 
with X appearing only at the next order. In addition, these correlators are independent of the 
unknown singlet constants mg and a at next-to-leading order in the quenched theory. It is 
therefore expected that these two-point functions will yield particularly clean determinations 
of (the quenched) F. 
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Our paper is organized as follows. In Sec. [21 we set up our notation and reiterate a few 
facts about the e-expansion of chiral perturbation theory. In Secs. 01 and 01 we present our 
results for flavoured vector and axial current correlators in the quenched theory. Finally, we 
present in Sec. 01 the analogous results for the full theory, and we conclude in Sec. 01 


2. Basic setup 


Let us consider QCD in a toroidal volume V of average length scale L = V 1//4 . We assume 
that the volume is large with respect to the QCD scale, i.e. FL > 1. As in infinite vol¬ 
ume the lightest degrees of freedom are the would-be Goldstone bosons of chiral symmetry 
breaking. Their interactions can be described in terms of a chiral Lagrangian, which can be 
systematically expanded in powers of the pion momentum and mass over the cutoff of the 
effective theory, A ~ AirF. To leading order in the chiral expansion, 


F 2 

C-xPT = -^-Tr 


dpUdurf 


-Ljy 

2 


UMe id,N f + M^U ] e~ id,N f 


( 2 . 1 ) 


where U € SU (Nf), 6 is the vacuum angle, and, to this order, F and E equal the pseudoscalar 
decay constant and the chiral condensate, respectively. For simplicity, we consider the quark 
mass matrix M to be diagonal, M = diag(m,..., m). 

In an infinite volume both the momentum and the pion mass are taken to be of the same 
order. For a finite L smaller than the inverse of the pion mass, this conventional chiral 
perturbation theory breaks down when mT, V < 1, because the momentum zero-modes of the 
pion fields become non-perturbative T and require a resummation to all orders. Gasser and 
Leutwyler provided a procedure to rearrange the chiral expansion appropriately [TJ . The first 
step is to factorise the zero and non-zero modes by writing 


U = exp 



U 0 , 


( 2 . 2 ) 


where Uq is a constant (in space-time) matrix of SU (Nf), and the fields £(x) parametrize 
the non-zero-mode manifold. The power counting rules are then those of the so-called e- 
expansion: 


F~C>(1), 0 M ~0(e), L~C>( 1/e), f ~ 0(e), m ~ 0(e 4 ). (2.3) 


Thus, the quark mass goes like four powers of the momenta, instead of two. 

The difficulty of the e-expansion is simply that the integration over Uq has to be carried 
out exactly. Such integrations can however be performed, at least numerically, and results 
for meson correlators have been computed in full QCD at next-to-leading order in the e- 
expansion MW- In some cases, it might be interesting to consider these observables also in 
sectors of fixed gauge-field topology (particularly when small physical volumes are used). The 
fixing of topology in chiral perturbation theory amounts to an enlargement of the zero-mode 
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manifold from SU(IV/) to and to the addition of the term (detUo) 1 ' to the weight in 

the zero-mode integral EH, but then the integration can even be carried out analytically. 

Unfortunately, lattice simulations with very light quarks are still restricted to the quenched 
approximation. Therefore, analytical predictions will also have to be obtained using quenched 
chiral perturbation theory. 

Two recipes have been proposed for dealing with the quenched limit. One is the so- 
called supersymmetric quenched chiral perturbation theory na us], where the cancellation 
of internal fermion loops is achieved by adding to the theory ghost bosons of spin 1/2, one 
for each fermion field. A flavour singlet meson field, To, also needs to be added. The other is 
the so-called replica method nzi, in which only the flavour singlet field To (whose non-zero 
momentum modes then correspond to Tr£ in Eq. CT ) is added, so that the manifold is 
U (Nf). Computations are carried out by keeping track separately of the N v external flavours 
appearing in the operators and of the Nf propagating dynamical flavours, and the quenched 
predictions are obtained by taking the limit Nf 0, for a fixed N v . 

While it is believed that the two methods are equivalent in perturbation theory, the sit¬ 
uation is somewhat trickier in the e-regime, where a non-perturbative definition of the in¬ 
tegration over the zero-mode manifold is required. In particular, the replica method has so 
far provided results for the zero-mode integrals only in terms of series expansions in a mass 
parameter (see, however, the recent developments in USD- Therefore, replica computations 
are usually transformed into supersymmetric ones at this point. 

In the supersymmetric formulation, the zero-mode manifold is obtained by enlarging \J(Nf) 
to a graded group. The simplest choice is to take XJ(Nf\Nf) with Nf = N v , although in 
principle any other choice with a larger Nf > N v should be equivalent. To leading order in 
the momentum expansion the supersymmetric quenched chiral Lagrangian is 


F 2 

^QxPT = -r Str 


d tl ud tl u- 1 


mTi 


Str 


UgU + U^UF 1 


+ 


m o 

2 N, 


■H + 


w:(W 2 ' 


(2.4) 


where Str denotes the supertrace, To = ^Str[—iln(U)] and Ug = exp(iOI]\r v /N v ). Here, In v 
is the identity matrix in the fermion-fermion block of “physical” Goldstone bosons and zero 
otherwise. At this order there are again the two couplings E and F, as in Eq. dm but also 
a series of new couplings, mg, a ,..., associated with the flavour singlet field, which cannot be 
decoupled in the quenched limit da CE]. 

As discussed in mm, the e-expansion can be set up for this theory as well, at a fixed 
topological charge v. The integral over the zero-momentum modes requires, however, special 
care. To obtain convergent answers, the manifold is in fact not taken to be U (N v \N v ) but what 
in the mathematics literature is called the maximally symmetric Riemannian submanifold 
G1(W|A/;) m ■ The resulting quenched integrals are known explicitly for N v = i [2H 12 a, 
and some of them also for N v = 2 ESI . 

In ref. {7( the scalar and pseudoscalar two-point functions were computed in the e-expansion 
to next-to-leading order. Below we determine the analogous vector and axial vector corre- 
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lators in this expansion, at the same relative order. The reader can find further technical 
details on the method in [? . 


3. Quenched vector and axial current correlators 


To define the observables considered at the quark level, consider for definiteness a convention 
where the Euclidean QCD Lagrangian is Ce = '07 /t T ) // '0, with 7 ^ = 7 ^, { 7 ^, 7 ^} = 2 5^, and 
75 = 7o7i7273 • The vector and axial currents can then be chosen as 

V£ 0) = #(x)T^ b 7 /x ' ip(x) , A“(x) = 7^75 ^(x) , (3.1) 

where are traceless generators in the physical flavour space of the iVy valence quarks. 
We will work with the conventional normalization, 

Tr[T^J = (3.2) 

The corresponding currents in the chiral theory can be obtained by coupling covariantly the 
pion field to the external vector and axial-vector 0 “ fields, and then taking functional 
derivatives with respect to those sources. More precisely, the above conventions imply that 
the partial derivatives in the chiral theory are promoted to covariant ones as 


d^U- 

► + i 

«-<)T^C/-iC/T^(n“ + a“), 

(3.3) 

and the currents become 




III 

PT 

III 

H 

Oh 

* 

(d£ x PT\ 
V dvfi ) 

-£tv + a^u)] , 

D“=° 

(3.4) 


( dC x p T \ 
v <9a“ / 

-i^Tr[T^(-a M C/H- 1 + a M 17- 1 [/)\ 

(3.5) 


Inserting here Eq. ( 1 T 21 ) . expanding in £, carrying out the contractions, and taking the 
replica limit Nf —* 0, the two-point functions for V“(x),M“(x) are easily determined at 
next-to-leading order. The general result can be written in a compact form by introducing 
t± = Tpj v ± U 0 Tft v U q 1 . and denoting 0^~{x) = V“(x), 0^ + (x) = M“(x), whereby (omitting 
contact terms) 

( 07 ( 1 ) oy<o))„ = 

- ^r <TT ( Tr K t ”’ *1 + t/ 6" 1 )])t,c/ 0 / d 4 2 - x)9 t A(z - 0) . (3.6) 

Here (• • -)^ Uo denotes the zero-mode average, 

, vAt _ fu 0 (' ' ')(det Uq) v exp(ptRe TV [?7o]) 

f Uo (det UqY exp(/i Re Tr [C/q]) ’ 1 j 
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while A(x) is the meson propagator without zero-momentum modes, 


A ( x ) = T7E 


0 ip-x 




(3.8) 


Furthermore, fi = mSh, while g! is an order-e 2 corrected value thereof [7]; in practice, 
however, // is not needed here, since the leading order results turn out to be pure numbers. 
The first term in Eq. (TTT7T1) is the leading, the second the next-to-leading order contribution. 
In the supersymmetric formulation the expression is formally identical, after the substitution 
Tr —> Str and the corresponding change of the group manifold for Uq. 

Before we go on to discuss the zero-mode integrations, let us emphasize two important 
points. First, none of the higher-order couplings of the chiral Lagrangian (the Lj’s of Gasser 
and Leutwyler) contribute at this order. This is in contrast to the standard p-expansion 
of chiral perturbation theory where, for example, L 4 contributes to the axial correlator at 
next-to-leading order MM- This shows the advantage of the e-regime with respect to the p- 
regirne in extracting the low energy constants E and F. Second, as shown by a comparison of 
Eq. (IM.lil) with the unquenched results (in Sec. 0, quenching does not modify the space-time 
dependence of these correlators because, remarkably, all potential double-pole contributions 
to the propagators cancel at both leading and next-to-leading orders. In other words, the 
additional constants ttiq and a of the quenched theory do not appear at all. This is in contrast 
to the scalar and pseudoscalar correlators, which show more important quenching effects (7]. 
This implies that the extraction of E and particularly F , which appears at the first order, 
should be easier from the axial and vector currents. 

We now return to the zero-mode integrations in Eq. (EH). As explained above, the replica 
method does not yet provide a recipe to compute them, apart from series expansions. How¬ 
ever, since the perturbative rules for the two formulations agree, we can at this point switch to 
the supersymmetric formulation. For a diagonal mass matrix, the results for the zero-mode in¬ 
tegrations are invariant in the linear transformation Tfj —> gT^g^ 1 , Tjy —> gT^g~ 1 for any 
g £ G1(N V \N V ), and should thus be proportional to the singlet tensor 5 ab . For (V“ (x)V f 6 ( 0 )) 
and (*4 ,“(.t)* 4^(0)), it is then sufficient to consider 

li ee (Str[(4) 2 ])^ o , 1\ ee (Str [(^) 2 ( U 0 + U t q" 1 )])^ , (3.9) 


where the superscript a is not summed over. The superscripts in the X’s, on the other hand, 
refer to the order in the expansion at which these contributions enter, as mentioned above. 

It is believed that we can now take, without loss of generality, the generators Tft to be 
in the simplest group allowing to define flavour non-singlets, i.e. N v = 2. The first integrals 
in Eq. EH can then be evaluated by a series of tricks explained in ref. 0. Using the explicit 
integrations over Gl(2|2) obtained in jH3) it is easy to show that for the normalization of 
Eq. EH- 

X° = 0 , Z\ = 2. (3.10) 
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Figure 1: Two types of external flavour contractions of two quark bilinears: (a) “connected” 
and (b) “disconnected”. 


The striking result is that the vector correlator vanishes identically at leading order. 

The integrals at the next order are not as easy to obtain, because third derivatives of 
the partition function for N v = 2 are needed, and they have to our knowledge not been 
computed before. We have however performed large and small mass expansions in the replica 
method using the results of m , and found that Zf vanishes to leading order in both of these 
expansions. (This, in turn, would imply that Z\ = (2/N v )(Tr [Uq + Uq 1 ))^ Uq .) This can 
hardly be a coincidence: it suggests that also Zf is identically zero. It should be noted that 
in the unquenched theory the non-singlet vector correlator is certainly non-vanishing urn- 

In order to understand this feature of the quenched approximation better, we present in 
the next section a new way of obtaining these results. In fact, we will argue that the vector 
correlator vanishes in the quenched approximation to all orders in chiral perturbation theory, 
both in the p- and e-regimes. The same insight provides a useful trick to obtain certain 
integrals in Gl(N v \N v ) from those based on the much simpler Gl(l|l). 


4. Quenched zero-mode integrals and an all-orders argument 

The argument relies on a relation that can be derived at the quark level, where the effect 
of quenching is immediately obvious. We start by considering the two-point function of any 
flavour singlet quark bilinear, O 0 (x). In general there are two ways of contracting the external 
flavours leading to two types of diagrams, shown in Fig. 0 In full QCD with Nf flavours the 
first one, which we shall call “connected”, will behave in terms of IVj-counting as Nj+0(Nj ), 
while the second behaves as Nj + 0(Nj). The idea is now to isolate the quenched “connected” 
contribution of one external flavour (N v = 1) by appropriately normalizing the approach of 
the full singlet quantity to the replica limit Nf —> 0, 

<O°(z)0°(»)>SS£2i ‘ 1 = lim T(0°(x)0°(„)r'. (4.1) 

Nf—*0 iv f 

The contribution of the “disconnected” diagram vanishes in this limit since it is O(Nj). 

Now that the quenched connected contribution for one flavour (N v = 1) is isolated, it is 
straightforward to get any non-singlet (N v > 2) quenched correlator in the limit of degenerate 
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quark masses. Indeed, the non-singlet correlators only contain the connected diagram and this 
diagram gives the same contribution no matter whether the two lines carry the same flavour. 
In other words, the N v dependence of this diagram is a trivial overall factor, which can easily 
be provided: the non-singlet quenched correlator is obtained by multiplying Eq. EB with 
the factor in Eq. EH’ 


(O a (x)O b (y)) quenched = Tr 




rpb 

1 N v 


lim -L(O 0 (x)O°(y))™ 

iV/—>0 iV f 


(4.2) 


Up to this point we argued at the quark level. However, since the relation of Eq. EH must 
be valid at any space-time separation, it must hold also in the effective chiral theory. The 
only subtlety involved in computing the matrix element on the right hand side of Eq. EH) 
is that the chiral Lagrangian for full QCD must include the singlet pseudoscalar field, since 
the decoupling of the singlet does not commute with the limit Nf —> 0. 

Assuming that the relation of Eq. EH indeed holds in the effective theory, it is easy to show 
that the vector correlator vanishes to all orders in the quenched chiral theory 2 . The reason is 
that, if we were to replace Tjy v —> Tjy = In { in Eqs. (13.41) . (13.51) and write U = exp(i2£/F), 
then = 0, while A ^ = —2 F d^Tr^. In other words, the vector singlet (i.e. baryon number) 
current vanishes identically 3 , as does then also the limit of Eq. (IP1) . On the other hand, the 
axial singlet current couples to the flavour singlet field, which cannot be decoupled in the 
quenched limit, leading to a non-vanishing contribution for the axial current. 

Since the non-singlet vector current has the quantum numbers of the p meson, we note 
that the result mentioned is consistent with, and in some sense provides justification for, the 
common lore that the quenched p “does not decay”. 

Let us stress that according to the argument just presented, the zero-mode integrals in the 
case of degenerate quark masses can all be obtained from the simplest case N v = 1. In the 
supersymmetric formulation this implies a set of very non-trivial relations between integrals 
in G^IVyllV,,) for various N v . In particular, we have checked that one can reproduce in this 
way all the results for the non-singlet scalar and pseudoscalar correlators of :7J, and those in 
Eqs. (13.101) . without the explicit integrals over Gl(2|2) of |23j : the only non-trivial integrals 
needed are in Gl(l|l). 

Similarly, the quenched singlet correlator for any N v can be easily computed, using the 
singlet computed for N v = 1, together with the quenched connected contribution obtained as 
in Eq. EH- The difference of these two quantities is the quenched disconnected contribution. 
However, once the connected and disconnected contributions for the singlet with N v = 1 are 
isolated, the full result for any N v is simply N v X (connected) + N^x (disconnected). From the 
point of view of the effective theory, this implies again non-trivial relations between certain 
integrals in Gl(-/Vy|IVy) and Gl(l|l), or the corresponding replica limits of unitary integrals. 

2 Obviously there are non-vanishing shorter-distance contributions. 

3 This is true to all orders, not only at the level shown in Eqs. 13.41 . 13.51 . 
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For completeness, we illustrate this point in Appendix A 


Using the above observations we can easily compute the integrals T± in Eq. (13.91) . and thus 
simplify the expression in Eq. (E3. To display the result, it is convenient to consider the zero 
components of the currents, and to project onto zero spatial momenta by the usual Fourier 
integration. Using the relations (up to contact terms) 


<9qA(x — 0) 


doA(z — x)doA(z — 0) 



(4.3) 

(4.4) 


where T is the temporal extent of the box and h\(r) 
obtain (no sum over o): 
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we finally 


Cr nched (*o) = / (V 0 a (x) V 0 »)„ 

J X 

ct m ( xo ) = [ MSM -4;(o))„ 

J-X. 


0 , 


(4.5) 


F 2 \ 2mA, u (n)T 2 f x 0 \ 
T 1 \ T / 


(4.6) 


where the only non-trivial zero-mode integral is the same as appears in the quark condensate 
obtained with Gl(l|l) (22) : 


E 




+ U+l(/l)k(/-l(/i) 


V 

H—) 
h 


(4.7) 


where I U ,K U are Bessel functions. 

To illustrate the results of Eq. (14.61) . we show in Fig. [21 the axial correlator normalized to 
the tree level value, for a realistic lattice setup (of volume V = TL 3 ). Note that the curvature 
at zero quark mass limit is quite sizeable, for v > 0. 

One apparent paradox in Eqs. (USD , (Hoi is the following 4 . By performing a global chiral 
rotation, ^ > ?// = exp(i<f> a T ar y^)ilj, which should be a good symmetry in the chiral limit at 

fixed volume, one can relate the vector and axial correlators through 


£<V»C(0)> = - TT £ ■ f “ tC (W-4t(0)<5 c s), (4.8) 

a a iV ^ a,b,c 


where 5 C S = dS/d(j) c is the variation of the action. The second term on the right is clearly 
proportional to the bare quark mass, since this is the only source of chiral symmetry breaking, 
and one might then expect that it vanishes in the chiral limit, leading to equal vector and 
axial correlators. 

4 We thank M. Liischer for pointing out this paradox. 
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L=24a T = 32a a' 1 = 2GeV 



Figure 2: The expression inside the square brackets in Eq. U.6] ), for T = 32a, L = 24a, 
a -1 = 2 GeF, and two masses, m = 0 MeV (dashed) and m = 5 MeF (dotted), as well as a 
topological charge v = 1,2. T/ie solid line is just the tree-level result. For the parameters in 
the chiral Lagrangian (which to this order are scale independent) we have assumed F = 93 
MeV, E = (250 MeV) 3 . 


Clearly, Eqs. (EH) and EH contradict this expectation. To be explicit we have, in the 
limit of small quark masses, 




S (? + 2 ^)+^ 3 ), 

£ (l - 7 + log (2/V)) ii + 0(v 3 ), n = 0 


(4.9) 


which shows that the axial correlator in Eq. EH does not vanish in this limit. This is however 
not in contradiction with Eq. EH, because the second term on the right does not vanish 
either. The reason is that when averages are considered in sectors of fixed topology, there are 
chiral singularities in the three-point function on the right hand side of Eq. EH that cancel 
the explicit factor of m, giving a finite contribution in the chiral limit, which exactly matches 
the difference between the vector and axial correlators. While the situation for v ^ 0 is very 
similar in the full and quenched theories, and these singularities can be clearly identified as 
the topological zero modes, in the case of v = 0, the behaviours are different in the two cases: 
the vector and axial correlators are the same in the chiral limit in the full theory while they 
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differ in the quenched case. This indicates a more singular chiral behaviour of the quenched 
theory, which was already noticed in the case of the quark condensate PHD]. 


5. Results in the full theory 


For completeness and comparison we present in this section the results for the vector and 
axial correlators in the full theory at fixed topology. They can almost be read from (5J, if not 
because there are a few places where identities were used for the zero-mode integrals that are 
not valid when averages are taken in sectors of fixed topology. After integrating over space 
and using the relations 

N f 


J_ = (TV [(i“ y ([To + Uq)])^ Uo = -f(TV [(£)' - (t a _)1)^ 0 , 


Jl = (Tr[(t a + )\Uo + tf 0 )})^ Uo = ^(TV[(t“) 2 -(t“) 2 ])^ o + -^(ReTVC/ 0 )^ o) (5.2) 


N 




we obtain (no sum over a): 


J1 + 


c-(z„) ^ / x m*)v;( 0 »„ = -A | 

Cw ^ JjAZ(x)A‘ 0 m« = + 


N f 






+ 


Afi T 2 
N f F 2 V 


where pi , A:oo are numerical factors ESI IS], and 

Jl ee <TV[(t“ ) 2 ]>^; c/0 = _ 2 + <Tr[C/ 0 ] Tr[c4])^ 0 ) , 


J -- = <Tr[(t“) 2 ])^' n = 


v,Uo 


N f ~ 1 


(Nj - {Tr[U 0 ]Tr[U^ Uo ) . 


Here again / jl = rnT,V, while its order-e 2 corrected value is PP 

Nf ~ 1 Pi 


»=^{ 1 + N f F 2 V 1 ! 2 J 
The average left over in Eqs. (15.51) . (15.61) is given by 0 


(Tr[U 0 ]Tr[U^ Uo =N f 


W) 


+ Nr 




+ 


1 s„(/z 0 
n' s 


- Nf 


V ) 2 


(5.1) 


(5.3) 


hl ^)(ReTvU 0 rA , (5.4) 


(5.5) 

(5.6) 

(5.7) 

(5.8) 


where £„(//) = (E/Afy)(ReTV Uo)^ Uo is the unquenched chiral condensate at fixed topology 
(see for instance JH]). Note that the unquenched £„(//)/£ depends also on Nf, although 
this has not been displayed explicitly. 
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As we see from Eq. <53>, the vector charge is conserved in the full theory just as in the 
quenched one, but it has now a non-vanishing value. It can be easily checked, however, that 
at, say, u ^ 0 but [i —> 0, one gets = 0 + O(Nj), J+ = 2 + O(Nj), j}_ = 0 + O(Nf), 
j\ = (2/A r /)(Tr [Uo + C/q]} + 0(Nf ), in accordance with the quenched results in Sec. Q]for 
'I [) _, X { \_ , X\_ . , respectively. In the full theory the vector correlator has also been computed at 
next-to-leading order in the p-expansion M and, being proportional to Nf, again vanishes 
in the replica limit Nf —> 0. 

6. Conclusions 

We have presented the results for flavoured vector and axial vector two-point functions in 
quenched chiral perturbation theory, at next-to-leading order in the e-expansion, correspond¬ 
ing to a finite volume and the vicinity of the chiral limit. Analogous results for unquenched 
QCD were already presented in Eli 

We find that quenching has a striking effect on the vector correlator, which can be argued to 
vanish to all orders in the chiral expansion. The axial correlator on the other hand depends 
at the leading order on the low energy constant F. and at the next-to-leading order also 
on E. However, it does not depend to this order on any of the couplings of the quenched 
theory associated with the singlet pseudoscalar field, m g, a ,..., nor on the L,’s of Gasser and 
Leutwyler. The measurement of the volume and mass dependence of this correlator close to 
the chiral limit would therefore permit the extraction of (the quenched) E and F , with a 
minimal contamination from higher order effects. 
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Appendix A. Examples of quenched zero-mode integrals for N v > 1 

To illustrate what can be achieved with the method discussed in Sec. 01 we show here how it 
helps to determine certain integrals with N v > 1 from known ones with N v = 1. We do this 
for the flavour singlet scalar and pseudoscalar correlators, computed for N v = 1 in [7] . 

Let us define 


S°(x) = i’(x)I Nv 'ip(x) , P°(x) = ^(x)I Nv i^(x) . (A.l) 
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The correlators of S°(x), P°(x) have two parts: a constant and a space-time dependent 
contribution. We consider here just the former. As in [7|, we denote the constants for S' 0 , P° 
by Cg and Cp, respectively. For N v = 1 it was found that, to lowest order, 


C° s = ESUm) + S 2 + ^ 

/i 


c% = 


sV 


(A.2) 

(A.3) 


I 1 ^ 

where £j,(/r) is the chiral condensate. Using now the prescription of Sec. 0 (Ed- Ijd.Ul l. we 
get the “connected” part of Cg by a U (Nf) computation, 




N f ^o Nf l 4 


= lim -1- ( E 2 AU 
N f ^o Nf l 1 


K(») 


+ N f 


Ev(fj) Y 


= EE'„(/x) - (A.4) 


Compared with the N v = 1 result in Eq. (TOll . this gives the full expression, 


sV' 


^quenched = + W 2 f E 2 + —\ 


(A.5) 


Similarly, for the constant Cp we get the “connected” part from a U (Nf) computation, 

. — | —/i'tv rrr„ _ 7rtr\2' 

Af/^o JVj 1 4 




= lim 


iV/—>0 Nf 


S 2 JV f 


1 _ v 2 Nf 

H £ ^u 2 






(A.6) 


Comparing this with the iV,j = 1 result, the general expression for any _ZV„ is seen to read 


[C£] 


fO i quenched N v T,T, u (fj.) 1V 2 SV 


(A.7) 


t-t l l 

Eq. (IXTl) is in agreement with what one gets from a chiral Ward Identity at fixed topological 
charge for N v = 2, as explained in [Jj. The space-time dependent parts of the two singlet 
correlation functions can be constructed for arbitrary N v in a similar way. 
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